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Abstract 

We analyze diffractive deep inelastic scattering within perturbative 
QCD by studying lepton scattering on a heavy quark target. Simple ex- 
plicit expressions are derived in impact parameter space for the photon 
wave function and the scattering cross sections corresponding to single 
and double Coulomb gluon exchange. At limited momentum transfers to 
the target, the results agree with the general features of the "aligned jet 
model". The color-singlet exchange cross section receives a leading twist 
contribution only from the aligned jet region, where the transverse size 
of the photon wave function remains finite in the Bjorken scaling limit. 
In contrast to inclusive DIS, in diffractive events there is no leading twist 
contribution to <jl/&t from the lowest order (qq) photon Fock state, and 
the cross section for heavy quarks is power suppressed in the quark mass. 
There are also important contributions with large momentum transfer to 
the target, which corresponds to events having high transverse momentum 
production in both the projectile and target rapidity regions, separated by 
a rapidity gap. 
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1 Introduction 



The physics of rapidity gaps in deep inelastic lepton scattering (DIS) has been 
qualitatively described using the "aligned jet" model [|T|. This model forms the 
basis for several QCD studies of diffractive DIS, where various assumptions have 
been made concerning the interaction with the target (Pomeron exchange ||, 
BFKL ladders ||, soft color fields Q]). In this paper we wish to test and clarify 
the physics of the aligned jet model by choosing a particularly simple target, 
which allows a full perturbative calculation. We believe that such an explicit 
(and straightforward) QCD calculation can serve to illustrate several points of 
principle, and help distinguish between the physics of the photon projectile and 
the target "diffractive structure function" . 

In the aligned jet model, the virtual photon breaks up into an asymmetric 
quark-ant iquark pair. In the frame where the target is at rest, one of the two 
partons, say the quark, takes nearly all the photon energy u, and forms the 
"current jet". The antiquark has then a finite energy of 0(Aqcd/%), even in 
the scaling limit, Q 2 , v — > oo with x = XBj = Q 2 /2m N b / fixed. The antiquark 
also has low transverse momentum p± = O(Aqcd) with respect to the photon 
direction, and scatters softly on the target. (Conversely, in the frame in which 
the target has infinite momentum, the antiquark is seen as a quark emerging 
from the target with momentum fraction x, and the soft scattering probability 
corresponds to the structure function.) 

Since the antiquark momentum is proportional to 1/x, it reaches quite high 
values in the small x range accessible at Hera. Thus, in analogy to hadron scatter- 
ing, also the antiquark scattering is expected to develop a diffractive (color-singlet 
exchange) part, resulting in a "rapidity gap" between the antiquark and target. 

The soft target scattering cannot be calculated in perturbation theory. In in- 
clusive DIS, one relies on the QCD factorization theorem to separate the hard and 
soft momenta in the process. For diffractive DIS, however, there is no analogous 
proof of the separation of hard and soft subprocesses, so that the analysis is more 
model dependent and may provide new insights into the structure of hadrons and 
nuclei (see for example Ref. 0). 

Here we wish to study a very simple model for the target scattering by assum- 
ing that the target is a heavy quark ||. This selects Coulomb gluon exchange 
as the dominant high energy interaction. We will consider both one and two- 
gluon exchange contributions in order to model the inclusive and diffractive DIS, 
respectively. Our approach can be viewed in two ways: 

(i) As a toy model, consistent with PQCD (in fact it is PQCD, as the process 
we consider is theoretically conceivable). Many of our results will be similar 
to calculations involving more sophisticated assumptions about the target 
scattering, for example in terms of BFKL gluon ladders || . We feel that it is 
useful nevertheless to explicitly demonstrate those features which follow just 
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from the photon structure through the aligned jet mechanism, by making 
the simplest possible assumption about the target scattering. 

(ii) As the lowest order of a more complete calculation which also takes into 
account gluon radiation and scattering. As noted above, in the target rest 
frame the antiquark is seen as part of the photon wave function; similarly, 
the (high energy part of the) gluon ladder to the target is built from gluons 
radiated from the initial quark pair. It is thus not unreasonable to expect 
the nature (and mass) of the target itself to be rather immaterial once 
sufficiently many gluons are considered in the photon wave function. 

At small x and, consequently, large longitudinal momenta p z ~ 0(Aqcd/%), 
the parton interaction times in the target (which are of (9(1 fm)) are short com- 
pared to their (Lorentz-dilated) lifetimes. Hence the scattering amplitude fac- 
torizes into a Fock state formation amplitude (the photon wave function) and 
a scattering amplitude for each Fock state [0, |8|, [|. Moreover, the scatterer 
can move only a negligible transverse distance ~ Lp±/p z in a soft scattering on 
a target of length L. Hence the scattering amplitude is approximately diago- 
nal in impact parameter space. We show explicitly how these features arise in 
the PQCD amplitude for one and two gluon exchange at leading order in 1/x. 
These features, which follow from general quantum mechanical principles, must 
be present to arbitrary orders in perturbation theory, providing an important 
conceptual and technical simplification of higher order calculations. 

Since hadronization is an inherently non-perturbative phenomenon, it is not 
possible to prove that the rapidity gap between the target and diffractive scatterer 
is maintained in the evolution to the hadronic final state. Nevertheless, one can 
take a step in this direction by calculating the distribution of soft gluon radiation. 
This is analogous to the PQCD "proof" of the string effect in e + e~ — > qqg — > 



hadrons JTT||. From the structure of the photon wave function we verify that 
soft gluon radiation is indeed dominantly emitted in the rapidity region between 
the quark and the antiquark. We expect this feature to survive a scattering with 
color-singlet exchange, leaving a rapidity gap between hadrons in the photon and 
target fragmentation regions. 

An immediate consequence of our analysis is that the dependence of the 
single-gluon exchange cross section on the photon energy v dominates the multi- 
gluon exchange contribution by a factor of 0(\ogv), a feature that appears also 
in QED (T^j. This logarithm is due to scattering at large impact parameters, 
which for a color-singlet target will saturate at impact parameters of the order 
of the target radius. 

Our model also shows that there are important contributions to hard diffrac- 
tive scattering from small impact parameters of 0(1/Q). This should show up 
in ep scattering as events where both the projectile and the target fragment into 
large transverse momentum particles, with a rapidity gap between the projectile 
and target fragments. 
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The purpose of this paper is to present a self-contained perturbative model of 
diffractive DIS which bypasses the uncertainties due to Pomeron physics. With 
this in mind, we shall be rather explicit and assume no methods beyond the 
standard arsenal of perturbative QCD. 

2 Scattering amplitudes on a massive target 

In this section we will state the kinematic properties of our model, and present 
the calculation of the relevant amplitudes. We shall be concerned with the QCD 
amplitude for 

j*(q) +t(m t ) -> q(pi) + q(p 2 ) + t(m t ) , (1) 

where we take the target t to be a heavy quark of mass m t . We shall calculate 
the process ([!]) in the small-x limitQ, at leading order in PQCD, and study both 
inclusive and diffractive DIS, the latter being defined by the requirement that 
the qq pair be in a color-singlet state after the scattering. In both cases, we 
demonstrate how the amplitude factorizes into the photon wave function (as 
derived in the Appendix) and the qq scattering amplitude. 

The calculation is greatly simplified by assuming a large target mass mt, which 
implies Coulomb gluon exchange and negligible energy transfer. The large target 
mass limit is thus a natural laboratory for perturbative studies of DIS and for the 
emergence of rapidity gaps. Since, as we shall see, much of the physics is encoded 
in the photon wave function, there is reason to hope that many features of the 
process will be target-independent and generally applicable to the physical case 
of a color-singlet, finite mass target. 

The space-time picture we have in mind, depicted in Fig. 1 for double- 
gluon exchange, is thus the following. A high energy photon of energy v and 
virtuality Q 2 fluctuates into a qq pair a long time (0(u/Q 2 )) before reaching 
the target. The pair subsequently interacts with the massive target via single or 
multiple Coulomb gluon exchange, with a limited total momentum transfer K. 
For a diffractive process we require the pair to emerge from the scattering in a 
color-singlet state, so that a rapidity gap between the pair and the target can be 
formed. 

We will work in the target rest frame in which the virtual photon is moving 
along the positive z axis, q = [y, 0j_, \/v 2 + Q 2 ). The small value of Bjorken x 
implies that the photon is off-shell by an amount which is small compared with 
its energy, 

q z ~v + Q 2 /2v . (2) 

1 We use the standard definition of x = Q 2 /2mNV, tun being the nucleon mass. All our 
results will depend only on the combination xtun — Q 2 /2v and are insensitive to the value of 
the (large) target mass m t . 
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Figure 1: A diagram contributing to the process e + t(m t ) —> q{pi) + q{P2)+t{m t ) 
with double-gluon exchange. 

We take all transverse momenta to be much smaller than the longitudinal ones. 
The 3-momenta of the quarks are parametrized as 

Pi = (Pi±,zq z ) , 

P2 = (P2±, (1 - z)<u + K z ) , (3) 

where = Kj_ — p±±. Energy transfer to the target is suppressed by powers of 
the target mass, so that energy conservation reads 

Pl± + ^ Pl± , A N 

Izv 2(1 — z)v 

This implies that the total longitudinal momentum K z transferred to the target 
is small, of order Q 2 jv or p\jv. 

With these kinematic preliminaries in mind, we turn to the evaluation of the 
scattering amplitudes. In particular, using the standard decomposition of the 
total amplitude given in the Appendix, Eq. (|60D, we concentrate on the amplitude 
for process (|l]) with a virtual photon of definite polarization, 

A(A;Ai,A 2 ) = e(A).j h (Ai,A 2 ) , (5) 

where Ai and A2 are the helicities of the quark and the antiquark, respectively. 

We first calculate the single Coulomb gluon exchange amplitude A\, which in 
many respects is analogous to the two-gluon (color-singlet exchange) amplitude 
A 2 . In particular, both amplitudes will turn out to factorize in impact parameter 
space into a product of the photon wave function (c/ Appendix) and the qq 
scattering amplitude. 
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2.1 One— gluon exchange 

The scattering amplitude is the sum of the two diagrams given in Fig. 2. Consider 
for example the first diagram, in which the Coulomb gluon is attached to the 
antiquark. In the limit of large mass mt the target contributes only a factor 2mt 
times the color factor (T a )cD- The fermion propagator in the upper vertex may 
be split into forward and backward propagating parts using the identity 



+ m 1 



p2 _ m 2 _|_ j e 2E 



E'-f — p ■ 7 + m Ej + p • 7 — m 
p° - E + ie p° + E — ie 



(6) 



where E = \/p + m 2 - This procedure splits the covariant diagram into different 
time orderings of the vertices, reproducing the results of time-ordered and light- 
cone perturbation theory. In the DIS limit, the fermion energies are > 0{{p\ + 
m 2 )2v/Q 2 ), hence for small x we may neglect the backward scattering term in 
Eq. (|6]), and the denominator simplifies to 

2 2 

2E(p — E) ~ _ak±£_ (7) 

z 

where 

e 2 = m 2 + z(l - z)Q 2 . (8) 
The amplitude for Fig. 2a is then 

A?(p 1± ,K ± ) = ee q g 2 (T") AB (T a ) CD K2 ^ + £2) ( 9 ) 

X J2 [«Ai (Pl)e(A) • -yVp(p 2 - K) Vp{p 2 - K)70^A 2 (P2)] • 

P 

Using the explicit expressions for the spinors given in the Appendix, we have in 
the high energy limit 

v p (p 2 - K)-f vx 2 (p2) = 2(1 - z)v5 PM . (10) 
The second diagram is treated similarly, giving 

AfW.KJ = ee q g 2 (T°) AB (T a ) CD l\ 

^ IPi_l ' 6 ) 

x u A i(pi)e(A) -~fv\ 2 (p2-K) , (11) 

A (2), v v 2 /™X /ryn \ 4:171^(1 ~z) 

A\'(p 1± ,K ± ) = - ee q g (T ) AB {T a ) C D 



K^p 2 ± + e 2 ) 
x u Xl (p! - K)e(X) ■ ~fv X2 (p2) ■ 
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Figure 2: Single-gluon exchange diagrams contributing to inclusive deep inelastic 
scattering on a heavy quark target. 
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The close relationship between these amplitudes and the photon wave function 
given in the Appendix, Eq. ( [T3f) , is already apparent. Fourier transforming to 
transverse distance space r_i_ between the quarks, by defining 

A(r ± , K ± ) = J ^±A( P± , K ± ) exp(ir ± • p 1± ) , (12) 

we find consequently that the one-gluon exchange amplitude 

Ai(r ± ,K ± ) = ig 2 2m t (T a ) AB {T a ) GD W 1 {v^ K ± )V^ M (z, rj , (13) 

is proportional to the photon wave function V^ 1 ' 2 (z, rjj of Eq. (|82|) . The eikonal 
factor 

H r l(r±|K±) = iz^J^ (14 ) 

is independent of the photon helicity A and of the quark helicities Ai, A2. Since all 
dependence on the momentum transfer Ki to the target is in W\, we may make 
a further transformation to target impact parameter (RjJ space by defining 

Wi(r±,R±) = |^|Wi(r ± ,Ki)exp(iR ± -K ± ) 
= ^lim[K (\K ± \e)-K (\K ± + r ± \e)} 
= ^logf^^l (15) 



The physical interpretation of our result, Eq. (0), is straightforward. The qq 
pair, which is in a color-singlet state before the scattering, forms a color dipole 
of size r_i_. This dipole cannot be detected by a gluon with |Kj_| < l/|rj_|, in 
which case the two diagrams tend to cancel. Thus in Eq. (|T5|) , W\ oc |r_|_|/|Rj_| 
when this ratio is small. In particular, the infinite Coulomb phases associated 
with the quark and the antiquark cancel, and the eikonal factor is their finite 
remainder. In QED this factor would simply exponentiate upon adding extra 
virtual Coulomb photons. In non-abelian QCD the situation is complicated by 
color dynamics. The simple structure of Eq. ( JIB| ) will, however, be preserved for 
more than one Coulomb gluon, in the color-singlet exchange channel. 



2.2 Two— gluon exchange 

Four diagrams contribute to the scattering amplitude via two-gluon exchange, 
corresponding to different attachments of the gluons to the quark and antiquark. 
The individual diagrams are infrared divergent, but their sum is finite because of 
the dipole cancellation discussed above. We regulate the infrared divergences by 
allowing only a finite time r to elapse between the two gluon exchanges. With 
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a more realistic model for the target, r would play the role of the target size, 
and we would expect r ~ 0(1/ Aqcd)- Note that the true time interval between 
the two gluon exchanges is shorter than either the lifetime of the qq fluctuation 
or its subsequent hadronization time |7|, [13| . Our infrared cutoff r implies, as 
one can see by translating to time-ordered perturbation theory, an extra factor 
1 — exp(— iK lz r) in the amplitudes, where K iz is the longitudinal momentum 
of one of the exchanged gluons. As we will see, the results are not sensitive to 
the specific value of the cutoff, which can be removed from the sum of the four 
diagrams. 

Let us consider the diagram in Fig. 3, where the two gluons with momenta 
K x and K 2 (K x + K 2 = K) interact with the antiquark line. Upon using the 
identity (^|) on the fermion line between the gluon exchanges we encounter the 
energy denominator 

AE = -K u ~ jz 1 , \e 2 + P 2 1± + z(K 2 1± - 2p 1± ■ K 1± )1 . (16) 

lz[l — Z)V L " ' " J 

Since K\ z > 0(1/ r) due to our cut-off, we have AE ~ — K\ z at high energy v. 




Figure 3: A two-gluon exchange diagram contributing to diffractive deep inelastic 
scattering on a heavy quark target. 

It is legitimate to take the target mass limit mt —>■ oo inside the loop integral, 
since the integral is ultraviolet convergent. Then the target contributes just an 
overall factor of 2m t and a color coefficient, as for one-gluon exchange. The 
upper vertex also simplifies in the high energy limit, and finally we find 

^(Pix.Kx) = ee^ 4 (T*T b ) AB (T b T a ) CD ^H^zA (17) 

Pi± + £ 
f <i 2 K 

x u Xl (pi)e(X) ■ jv X2 (q - Pi) J J(K 1± , K 2 ±, t) , 



S 



where J is the integral over the longitudinal component of the loop momentum. 
Using the fact that K z = K\ z + K 2z is negligible, the longitudinal momentum 
integral can be computed exactly, 

r +oo dPC^ 1 - fi-^ T 

J (a, b, t) 



2tt (Kl + a 2 )(Kl + V)(K lz 



2 b 2 



-br 



b 2 



Notice that for large values of r this becomes simply 



it 



18) 



J (a, b, t) 



(19) 



2a 2 b 2 
independent of r. 

The other three diagrams can be treated in the same way. The longitudinal 
momentum dependence is encoded in the same integral, Eq. (18), for all of them. 
We find 

,(1,2), v v 4 frpbrpa^ frprp, fCPK 1± %Tl f V z{\ ~ £) 

A\ ; (pi±,K_l) = - ee q g (T T ) AB (T b T a ) CD 



(2tt) 2 (p 1± -K 2± y+e< 
x u Al (pi - K 2 )e(X) ■ -yv\ 2 {p2 - K 1 )J(K 1± ,K 2 ±,t) , 



4 2ll) (Pi±,K ± ) 



4 2 ' 2) (Pi±,K ± ) 



ee q g 4 (T a T b ) AB (T b T a ) CD J 



d 2 K 



i± 



(27T) 



lPl± 



K l± ) 2 + e 2 



x u Al (pi - Ki)e(A) • jvx 2 (p2 - K 2 )J(K 1± ,K 2± ,t) 

Arritvzil — z) 



ee^ 4 (T T a ) A B (T b T a ) CD - , 

P 2 ± + £ 



x u\i(q -P2)e(A) • 7^a 2 (P2) 



ci 2 K 



i± 



(27T) 5 



J(K 1± ,K 2± ,t) . (20) 



There are now two independent patterns of color flow, corresponding to singlet 
and octet exchange between the pair and the target. Since we are interested in 
the emergence of rapidity gaps, we project the amplitudes in Eqs. ( |17D and ( p0|) 
onto singlet exchange by summing over the colors of the target quark. This is 
crucial to the present argument since only upon performing this projection do 
the four diagrams acquire the same weights (as in QED), and the eikonal factors 
in impact parameter space are similarly reconstructed. Upon summation over 
target color the color factors of all diagrams in fact reduce to that of the first 
one, 



(T a T b ) AB (T b T a ) CD S 



CD 



-CfSab 



(21) 



where C F = (N 2 - 1)/2N. 

(1 2) (2 1) 

The transverse momentum integrals associated with A 2 ' and with A 2 in 
Eq. (|20"D are more complicated than the others, since one more propagator de- 
pends on the loop momentum. This difficulty is however bypassed by once again 
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turning to coordinate space, as in Eq. (|I2]) . After the target color summation of 
Eq. ([21]) we find that the complete two-gluon color-singlet exchange amplitude, 



A 2 (r ± ,K ± ) = i<? 4 2mAC F 8 AB W 2 (r ± ,K ± ) Vf> x *(z,r ± ) 



(22) 



is, just like the single-gluon amplitude of Eq. (|13"D , proportional to the photon 
wave function V^ 1 '* 2 ^, rj_) of Eq. Q82p. The two-gluon eikonal factor is 



W 2 (r ± , K ± ) = j ^±J(K 1X , K 2± , t) (l - e ir - K -) (l - e 1 



ir i Kn 



(27T) 



(23) 



where K 2 = K — Kj. The integral is now finite, due to the extra suppression 
for small K±± and K 2 ± provided by the eikonal numerators. It is thus possible to 
let the cutoff r — > oo, and use for J the simple expression given in Eq. (^). As a 
result we see that the two-gluon eikonal factor is a convolution of two one-gluon 
factors, Eq. (0), and can be written as 



W 2 (r ± ,K_ 



i f d 2 K, i 1 - e ir ^' KlJ 



K 2 



2 J (2tt) 2 
x (27r) 2 5 2 (Kj 



rf 2 K 2± 1 - e ir ^- K ^ 



(27T) 



/T 2 



K 1± -K 2± ) . 

Transforming M^2(r±,Kj_) to Rj_-space as in Eq. fll5]) gives 

W 2 (r ± ,Rj.) = ^[^(r ± ,R ± )] 2 

R± + r±f 

|R±I 



(24) 



8tt 2 




(25) 



The fact that the two-gluon eikonal factor W 2 is the square of the single 
gluon factor W\ is a consequence of the conservation of the impact parameters 
r ± ; R_i_ during the scattering, and will clearly generalize to an arbitrary number of 
exchanges in the singlet channel. It has an important consequence for the energy 
dependence of the cross section, which was first noticed in QED by Bethe and 
Maximon fll2"| . The cross section for single-gluon exchange involves an integral 
which is logarithmically divergent at large impact parameters R^ to the target, 



d 2 RjWM 2 oc 



^/Q 2 d 2 R i 



R 



(26) 



This results in a cross section which grows logarithmically with the projectile 
energy. For two or more gluon exchanges the integral is on the other hand 
convergent at large Rj_, and no logarithm is generated. 
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3 Cross— sections 

Our results for the one and two-gluon exchange amplitudes 7*t — > qqt (||) can be 
summarized (see Eqs. ( |T3"D and (|2"2])) as 

A n {r ± , R ± ) = c n m t W n {r ± , R ± ) V{z, r ± ) , (27) 

where r_i_, R^ are the transverse distances between the produced quarks and 
between them and the target, respectively. The photon wave function V(z, t±) is 
given in Eq. (|82|) and the eikonal factors W n (r±, Rj_) for n = 1, 2 gluon exchange 
in Eqs. ( |15D and (|25|), respectively. The coefficients c n , 

ci = 2ig 2 {T a ) AB {T a ) C D , 

c 2 = ig A C F 5 AB , (28) 

reflect the color structure. 

The spin-averaged square of the full et — > eqqt amplitude T n can be expanded 
as in Eq. (|65|) , 



|T n | 2 = ^^||c n iy„(r ± ,R ± )| 2 z(l-z)F(^r ± ) , (29) 



with 



F(z,v ± ) = l[l + (l- y Y]{[l-2z(l-z)]e 2 K 2 1 (er ± ) + m 2 K 2 (er ± )} 
+ 4(l-y)[z(l-z)} 2 Q 2 K 2 (er ± ) 

- 2(1 - y) cos(2^)2(l - z)e 2 K 2 (er ± ) , (30) 

where we have used the explicit expression for the photon wave function derived 
in the Appendix (see also Ref. ||). Here ip is the azimuthal angle between the 
lepton plane and the interquark separation rj_. 
The cross-section is given by 

da 1 4m w r dz 



Q WS^ = 1 w^> dK±dpi± |r|2 ' (31) 

where m t is the target mass. We may thus define a "cross-section" in impact 
parameter space as 

Q A ^, d 7^ = ^A—, t^z [ d 2 r ± \c n W n (r ±J R ± )\ 2 F(z,T ± ), (32) 



dQ 2 dxd(pd 2 R ± (2vr) 4 xy 2 Jo 

whose integral over R^ gives the usual DIS cross section. In the Bjorken scaling 
limit ( |63"D this cross-section should (at leading twist) be independent of Q 2 . 

In considering the properties of Eq. (|32"D it is important to distinguish the 
regions of large and small impact parameters R±, corresponding to small and 
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large momentum transfer to the target, respectively. Most of the HERA data on 
rapidity gaps is selected for small momentum transfer, with the target (proton) 
diffracting into a low mass system |14|. Large momentum transfers implies par- 



ticles or jets with large transverse momentum also at target rapidities. There is 
experimental evidence for this type of diffractive events as well [p 



3.1 Low momentum transfer to the target: R± ^> 1/m 

The Bessel functions Ko t i(er±) in Eq. (|30| ) limit the size of the quark pair to 
r\ < min{l/m 2 , l/z(l — z)Q 2 }. Hence for R± ^> 1/m we have also R± 3> r±, 
and the eikonal factors ([H)| ) and fl2l>D, integrated over the azimuthal angle ipn 
between Rj_ and r_|_, are approximately 

2n / r ^\2n 



d<p R \W n \ 2 = a n[-j^) , (33) 

where a\ = l/An and a 2 = 37r/(47r) 4 . The integrals over in Eq. ( |3"2"| ) can then 
be done using 

' drir?K](er ± ) = ^ , (34) 

where the bj n (j = 0, 1; n = 1, 2) are numerical constants. Altogether the cross 
section becomes 

Q 



i da n 



dQ 2 dxd(fdR 



2 



a e q Q 2 | 2 a 

R x >l/m 



4(2tt) 3 x?/ 2 1 n| Rf 



x ^ {l[1+(1 _ y)2] [ (1 _ 2 , (1 _, ))£261 „ +m260 „ 



+ 4(l-y)[z(l-z)} 2 Q% n -2(l-y)cos(2<p)z(l-zyb ln } . (35) 



There are two regions of the final z-integral in Eq. (p5|) that should be con- 
sidered separately. 

(a) z = 0{m 2 /Q 2 ), r± = 0{l/m). 

This is the "aligned jet region" ]I|], relevant to the lowest order DIS process 
7*g — > q. In this region e is finite (0(m)), so the exponential suppression 
in the Bessel functions in Eq. ([34]) forces r± to be of order 1/m. Thus in 
the scaling limit one of the quarks has a finite momentum, zv ~ m 2 /2m^x, 
and the transverse size of the quark pair also remains finite. Only terms in 
the integrand of Eq. ( |35"D that are of O(z ) contribute. Hence the last two 
terms, which arise from photons with A = and from A = ±1 interference, 
respectively, can be ignored, and the leading twist cross section is isotropic 
in azimuth. Since e is finite, this region contributes both to the single and 
to the double-gluon exchange cross sections {n = 1,2). This is as expected, 
since attaching several gluons to the quark pair does not change the scaling 
behavior of the cross section when the transverse size of the pair is fixed. 
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(b) z = 0{\), r ± = 0(l/Q). 

In this region both quarks take a finite fraction of the photon momentum, 
so that e = 0(Q), which in turn forces the transverse size of the pair to 
decrease with Q. The photon fragments thus typically have large transverse 
momentum. Since e oc Q, the contribution from this region to double- 
gluon color-singlet exchange (n = 2) is suppressed by a factor l/Q 2 , ie, 
it is of higher twist. This is due to the compactness of the quark pair, 
which suppresses gluon attachment. In the case of single gluon exchange 
the suppression is compensated by the larger phase space available in z, 
compared to case (a). In the infinite momentum frame, the contribution of 
region (b) with n—1 corresponds to the process j*g — > qq in inclusive DIS, 
which is suppressed only logarithmically by the running coupling a s (Q 2 ). 

We conclude that, for moderate momentum transfer to the target, the color- 
singlet (n = 2) cross section fl3"o] ) gets a scaling contribution only from the aligned 
jet region (a). After the final z integration, we find 



Q 4 <h " 



dQ 2 dxd(fdR A 



* 2 e 2 q \c 2 \ 2 
4(2tt) 3 xy 2 

^^^[1 + (1 " ?/) 2 ] (2&i2 + &02) • (36) 



Notice that, according to Eq. (|3q), the heavy quark diffractive cross-section de- 
creases with the quark mass m as 1/m 2 . This is analogous to the suppression 
of diffractive production in the symmetric region (b) above. The size of a heavy 
quark pair is small, r± = 0(l/m), suppressing the coupling of multiple gluons. 
However, heavy quarks can be produced at leading power in the quark mass 
through higher Fock states of the photon, such as qqg. The two-gluon exchange 
can then occur off the soft gluon in the Fock state, if it has a finite momentum 
and thus also a finite transverse size distribution in the scaling limit (see Section 

It is also interesting to compare our analysis with the perturbative QCD 



calculations |T^] for diffractive vector meson leptoproduction: 7*p — * Vp. In that 
case the leading contribution to the cross section arises from longitudinally- 
polarized photons and the symmetric regime z = 0(|), r± = 0(1/ Q). Because 
of the color cancellations the exclusive diffractive cross section is suppressed by 
a factor l/Q 4 relative to the Bj or ken-scaling rapidity gap rate. 



3.2 Sum over target momentum transfers 

In the remainder of this section we will focus on color-singlet exchange. Since 
the eikonal factor W^rj^R^) according to Eq. (|25| ) depends only on the ratio 
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|R_i_ + rj_|/|Rj_|, an integral over all impact pararameters Rj_ in the cross section 
l2|) can be readily done. By dimensional analysis 



J d 2 R ± \W 2 (v ± ,R ± )\ 2 = d 2 r 2 ± 



(37) 



where d 2 is a calculable numerical constant. When unconstrained, the typical 
value of R± in the cross section (|32|) is always of 0{r±). Note that the power of 
r± in Eq. fl37D differs from that of Eq. (|33|) (valid for R± 3> r±). In terms of the 
constants bj n of Eq. ( |34"D we thus find the Rj_ -integrated cross section to be 



Q 4 



do-y 



7ia 2 e 2 Q 



dQ 2 dxd(f 



(2vr) 4 xy 2 



l2 , /"§ dz 
c 2 | d 2 / —7 







x {§[1 + (1 - yf] [(1 - 22(1 - z))e 2 b n + m\ 
+ 4(1 - y)[z(l - z)] 2 Q 2 b m - 2(1 - y) cos(2<^)z(l - z)e 2 b ll ) 



(38) 



In region (a) of the z integral in Eq. (|38D (z = 0(m 2 /Q 2 ), cf Section |3.1| ), r± 
(and hence also R±) is finite in the scaling limit. The scattering then dominantly 
proceeds via moderate momentum transfer 0(l/R±) to the target. On the other 
hand, for finite values of z in Eq. (|38|) the typical values of R± are of 0(1/Q), 
implying large momentum transfers to the target. The corresponding inclusive 
DIS process is then j*t — > qqt, where for a realistic target t would be a light quark 
(or gluon) constituent. Now the two-gluon exchange process is not suppressed 
since the typical hardness of the gluons is commensurate with the size of the 
qq dipole. This hard diffractive process (like its inclusive equivalent) is only 
suppressed by powers of the running coupling a s (Q 2 ). 

The z integration in Eq. ( |3"8| ) can readily be performed, and the leading twist 
contribution extracted. We find 



Q 4 



da 2 



2^2 



na e 



l C 2| 



dQ 2 dxd(f 



(2tt) 4 xy 2 



: [l + (l-2/) 2 ] 



An 



1 - lot 



Q 2 



+ (l-i/)(26oi-6ucos(2^)) 



(39) 



Notice that for hard diffractive events the azimuthal distribution is not isotropic, 
and a logarithmic dependence on Q 2 is generated, since a much larger range in 
p± is now available. 



3.3 Discussion 

To summarize, there are two quite distinct kinematical regions contributing to 
the leading-twist color-singlet exchange cross section, Eq. (^). For small mo- 
mentum transfers to the target, or equivalently for large impact parameter R±, 
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our calculation provides a concrete realization of the aligned jet model PJ: asym- 
metric quark pairs with finite transverse size dominate in the scaling limit, with 
the scattering occurring off the slow quark. In this case the cross section is purely 
transverse and independent of azimuth. If, however, we allow for large momentum 
trasfers to the target, we find that the color-singlet cross section receives lead- 
ing twist contributions from symmetric and narrow pairs. These contributions 
have a distinct azimuthal dependence and rise logarithmically with Q 2 . They will 
appear as events with particles (or jets) of large transverse momentum in both 
the current and target fragmentation regions, which balance each other and are 
separated by a large rapidity gap. There is already experimental evidence for 
such events from ZEUS ||15|| . 

Our diffractive cross section fl39|) has no contribution corresponding to events 



with large transverse momentum jets only in the current fragmentation region, 
with low momentum transfer to the target. Such contributions should occur at 
higher orders in QCD, for qqg Fock states of the photon (see Section 4). 

As we have seen, the low-x cross section is particularly simple when expressed 
in impact parameter space. Impact parameters are not directly measurable, since 
events are observed in momentum space. For integrated cross sections this is not 
an issue, as an integral over all impact parameters is equivalent to an integral 
over all transverse momenta. On the other hand, more differential predictions 
concerning, eg, the dependence of the cross section on the mass M 2 of the diffrac- 
tive system, and on the size of the rapidity gap At], are less straightforward. At 
this stage we only make the following qualitative remarks. 

The dependence of the diffractive cross section fl32|) on the mass M of the 
(qq) system is usually parametrized through the variable^] 

O 2 e 2 - m 2 
3 = ^ = - - (40) 

where e 2 = m 2 + z(l — z)Q 2 , as before. The rapidity gap r] extends from the 
target (which is at rest and thus has zero rapidity) to the slow quark, and (for 
z < 1/2) is given by 

(1zv\ , , 

At? = log , (41) 

Vrnii/ 



where mx± = \Jp\x_ + rn 2 . Notice that we use a frame where the 7* is moving 
along the z-axis. 

We can estimate the average size of the rapidity gap in the various kine- 
matical regions contributing to the cross section from the uncertainty relation 
corresponding to Eq. (0), 

Pl± ~— . (42) 



2 In the framework of Pomeron exchange models, j3 is the fraction of Pomeron momentum 
carried by the struck quark. 
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In the aligned jet region we have (as discussed above) (r±) ~ 1/m, and hence 
(Pii) ~ m - Since in this region z ~ m 2 /Q 2 the typical rapidity gap will be 

(nrn \ 
• (43) 

For large momentum transfer to the target, on the other hand, region (b) con- 
tributes to the cross section, so that (r±) ~ 1/Q, and consequently (pi,j_) ~ Q. 
Furthermore, in this region z ~ 1/2, so that the size of the gap will now depend 
on Q. We get 

A77 ~ log . (44) 

\mNXj 

The typical value of (3 remains finite in both regions, as it should. However we 
expect region (b) to give a relatively larger contribution at small values of (5, since 
the diffractive system contains large p± jets. Furthermore it should be kept in 
mind that higher order contributions from gluon radiation will tend to increase 
the cross section at low (5. 

The estimates (|43|), (fH|) for the size of the rapidity gap are at the parton level. 
Hadronization may to some extent fill in the gap. In the next section we shall 
address this question by studying the angular distribution of soft gluon radiation. 



4 Soft Gluon Radiation 



In the case of color-singlet exchange one expects that hadrons will be produced 
mainly in the rapidity interval between the quarks, leaving a gap between the 
slower quark and the target (here at zero rapidity). While this hadronization 
pattern cannot be demonstrated using perturbative methods, it may be corrobo- 
rated by studying the distribution of soft (but perturbative) gluons [fTOfl . This is 
analogous to the successful description ]TTJ of the "string effect" in e + e~ annihi- 
lation. 

Just as in the qq case above, we expect the gluon distribution in the final 
state to be largely determined by its distribution in the photon wave function. 
At lowest order, the qqg component of the photon is given by the two diagrams 
in Fig. 4. Using the prescription (fflD and the identity (j6|) to select the dominant 
time ordering, we find 



V(qqg) 



I . - u(pi) 



AE 



- e(g)-7 



1 £ fc i7° - (k + pi) • 7 + m 



1 £ fe2 7 C 



2-Efei E, 
(k + p 2 ) • 



Ei 
- m 



2E, 



A- 2 



E~ — Ei — E t 



k2 



E k i + ie 



v(p 2 



-7 



(45) 



where 

Ey El - 



E 2 - Ikl 



+ ie. 



Pi + m 2 , 



E, 



hi 



(k + pj) 2 + m 2 and AE 
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In the soft gluon limit, k z <C pi z , p 2z and |kj_] <C |pi_i_|, |p2±|, the qqg Fock 
state amplitude (|45|) factorizes into a product of the qq amplitude and a gluon 
emission factor, 



V{qqg) = V{qq)gT% B -^s\k) 



El _ El 

E\ E2 



(46) 



with V(qq) given by Eq. (|73|). The gluon angular distribution is thus governed 
by the factor 



P(k) 



— E 

(AE) 2 t 



El _ Ei 

Ei E2 



(47) 



where the sum is over the two transverse polarizations available to the gluon. 

To evaluate Eq. (|47|), let us for clarity choose the frame so that the gluon 
is emitted in, say, the x — z plane, ie, = (k,ksm8,0,kcos8). Within the 
approximations discussed in Section 2, we find that the sum over polarizations 
yields 

2 cos 2 d{plf + {p\f 



E 

A 



El _ El 

Ei E2 



v 2 z 2 {\ — z) 



while the energy denominator is given by 



» ^ pi + m , .„ 

- AE = mArX + r + fc(l - COS I 

2zA2ll — Zl 



(4* 



(49) 



For a given qq configuration (p±,z), and a given (soft) gluon energy k, the 
angular distribution of the gluon will be given by 



P{k) oc 



where 



A 



cos 2 0{pi) 2 + (p 
[1 + A(l - cos 1 

2vz(l - z) 



y \2 



l2 ' 



Pi + e 2 



-k 



In terms of the gluon rapidity 



1 . 1 + cos 9 

- log 

2 s l-cos# 



log tan 



e 



and of the (slow) quark rapidity r\ q given by Eq. (p]) we have 

2A; e^ 



1 + Ail - cos 1 



mj_ 1 + e 2r if> 



(50) 



(51) 



(52) 



(53) 



This denominator will suppress the probability (|50| ) of gluon emission unless the 
gluon rapidity is sufficiently high, 



27] ^ l ° g +Vq 



(54) 
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The softest gluons with k ~ m± thus have rapidities rj g £ ^ 3 /2. Harder gluons 
with k ~ zv will be emitted at rapidities between the fast and the slow quark, 
T} g £ i] q , outside the rapidity gap to the target. 

Thus the gluons associated with the photon Fock state tend to have rapidities 
in the photon fragmentation region, since this minimizes the off-shellness of the 
photon wavefunction. 

The above argument concerns gluons that are "preformed" in the photon wave 
function before the target interaction. The two (Coulomb) gluon interactions 
with the target occur during a very short time interval from the point of view 
of the photon. Hence, although the projectile system is momentarily a color 
octet after the first gluon exchange, there is no formation time available for (fast) 
gluon emission before the second gluon is exchanged. Thus we do not expect that 
gluons emitted in the short interaction interval can fill more than a limited part 
of the rapidity gap, close to target rapidities. 

We have in this paper neglected higher order contributions where the diffrac- 
tive scattering occurs off a slow gluon, rather than a quark, in the photon wave 
function. The amplitude will still factorize as in Eq. fl27|) into a Fock state ampli- 
tude V and a gluon diffractive scattering amplitude W. At lowest order the gluon 
Fock state amplitude should approximately be given by V(qqg) of Eq. (f46|). The 
quark configuration will now be a typical one, with z = (9(1/2) and p\ = 0(Q 2 ). 
The gluon must in this case contribute significantly to AE of Eq. (0^), having 
a finite k± in the scaling limit to maintain a constant diffractive (two-gluon ex- 
change) cross section. Just as in our analysis of quark diffractive scattering, this 
implies k z = 0{k\/m N x). The probability for such gluon Fock states scales like 
1/Q 2 , as can be seen by using Eqs. HSI and in the small 9 limit, 



P± [ A \h2 [ k ^ /mNX dk z 1 _ n( A \ 

as v 2 z\\-zfh dk± J 17 (m N x + k 2 ± /2k z ) 2 - U[as Qt } ■ (55j 

Since the gluon diffractive cross section will be independent of Q 2 due to the 
finite transverse size of the gluon distribution, we have a scaling contribution to 
hard diffractive scattering. By including Fock states with multiple gluons we may 
expect to gradually build up the gluon ladder or "hard Pomeron" j3|. 



5 Conclusions 

In this work we have shown a way to systematically analyze, fully within per- 
turbative QCD, the physical origin and dynamical dependence of rapidity gaps 
in deep inelastic lepton scattering. Our major simplifying assumption was the 
use of a heavy quark target so that the physics of the PQCD Pomeron could be 
directly identified with Coulomb gluon exchange between the heavy target and 
the constituents of the virtual photon. Color-octet exchange from one or more 
Coulomb gluons coupled to the target is expected to generate a radiative pattern 
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in the final state which occupies the entire rapidity interval between the target 
and virtual photon fragmentation regions. On the other hand, color-singlet ex- 
change arising from two or more Coulomb gluons corresponds to an absence of 
soft radiation in the central rapidity region. 

Our calculation shows explicitly the conservation of impact parameters in DIS 
at small x 0. This leads to a very simple factorized structure of the scattering 
amplitude, with the two-gluon eikonal factor W 2 (Eq. (^5|)) given by the square 
of the single gluon factor W\ (Eq. (|15D). 

Due to our use of perturbation theory, the explicit expressions we obtained for 
the scattering amplitudes incorporate general features required by field theory. 
It is, for example, a well-known feature of QED that single (Coulomb) photon 
exchange in Bethe-Heitler lepton pair photoproduction gives a cross section that 
increases logarithmically with photon energy v, 



The logarithm is due to an integration over large impact parameters R±, the 
upper limit being kinematically given by R± < l/K™ m = m 2 /2v. In our approach 
this logarithm is a direct consequence of the fact that W\ oc r±/R± for large 
target impact parameters R±, cf Eq. (|26|) . The fact that multiple Coulomb 
photon exchange does not give logarithmically enhanced contributions JT2| then 
follows directly from Eq. (]25f) , W 2 oc W 2 oc rj_/Rj_, and its generalizations. 
This difference in impact parameter dependence shows that the inclusive and 
diffractive DIS processes are dynamically distinct, and that their ratio cannot be 
characterized by a single color-dependent constant. 

In QED the logarithm of Eq. (|5"6|) saturates at R± = 0(1 A), where the 
scattering becomes coherent over the entire electrically neutral atom. In QCD 
the saturation for gluon exchange will occur at the color confinement radius 
R± = 0(1 fm), where perturbative methods cease to apply. 

In analyzing the structure of our diffractive cross section (^) we found it 
important to distinguish between the case of fixed (small) momentum transfer K 
to the target, and scattering involving a target momentum transfer that increases 
with the virtuality Q 2 of the photon. 

For small momentum transfers K we found that only the "aligned jet" region 
contributes to diffractive scattering in the scaling limit. This region is charac- 
terized by one of the quarks taking nearly all of the photon momentum, and 
by the transverse size r± of the quark pair remaining finite in the scaling limit. 
Hence multiple gluons can couple to the pair at leading twist. This is in contrast 
to the "symmetric jet" region where the vanishing size of the pair allows only 
a single gluon to couple. Since the size of a heavy quark pair is r_|_ < 0(l/m), 
two-gluon exhange is similarly suppressed in heavy quark production, leading to 
a diffractive cross section proportional to 1/m 2 . Analogously, the production of 
large pj_ jets at the photon vertex is suppressed by l/p± in diffractive scattering. 




£n — 

mi 



(56) 



20 



It will be important to search for these effects experimentally. These features are, 
however, true only at lowest order, when transverse gluon radiation is neglected. 
Most of the HERA data has been selected for low-mass target diffraction, 



which kinematically allows for a long rapidity gap between the photon and target 
fragments. It is clearly of interest to study also whether the hardness of the vir- 



tual photon can be transmitted to the target fragments in diffractive events [[15 
This would not be expected in a model where the scattering occurs via "soft" 
Pomeron exchange, similar to that observed in hadron scattering at low momen- 
tum transfers. In our approach we found, to the contrary, that the contribution 
from small target impact parameters R± is quite important and only suppressed 
by powers of the running coupling a s (Q 2 ). In inclusive DIS the corresponding 
subprocess would be classified as j*t — > qqt, where t can be any parton (light 
quark or gluon) of the target. Once = O(r^), the exchange of multiple glu- 
ons to the qq pair is not power suppressed no matter how small r± is, and the 
diffractive process can proceed at leading twist. In this large momentum transfer 
region perturbative calculations are actually the most reliable. 

We also considered the modifications brought to the above picture by radiative 
gluons. Soft gluons emitted in color-singlet exchange processes will mostly occupy 
the rapidity region between the quarks produced at the photon vertex, thus 
leaving intact a central rapidity gap extending to the target fragmentation region. 
At higher orders there are also contributions from photon Fock states where the 
quark pair is in a transversally compact "symmetric jet" configuration while a 
soft gluon maintains a finite transverse distance to the quark pair in the scaling 
limit ||. The soft gluon can scatter diffractively (through two-gluon exchange), 
allowing heavy quarks and large p± jets to be produced at the photon vertex 
at leading twist. These contributions tend to have large diffractive mass, and 
thus contribute to the region of low (3 = Q 2 /(Q 2 + M 2 ). They are, in fact, the 
beginning of the gluon ladder which presumably generates the BFKL Pomeron 
[0 . It should be possible to extend the present calculation in this direction - the 
more rungs in the ladder that are considered the less important should be our 
assumption of a heavy quark target. 
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Appendix 

Here we wish to collect some useful formulas concerning electroproduction and 
the virtual photon wave function. 



Al. Virtual Photon Factorization. 
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Figure 5: The amplitude for deep inelastic lepton scattering, £ + h £' + X. 



In the notation of Fig. Al, the lepton scattering amplitude is 

1 

Q 



T(£h^£'ti) = — j?(-g, u )f h , (57) 



where = e ■u (J '(£ / )7 /i 'U -(£) is the lepton current and j% is the target (hadron) 
vertex. We choose the target (h) rest frame where the virtual photon momentum 
is aligned with the z-axis, 



q=(u,0 ± , yV + Q 2 ) • (5* 
Gauge invariance, q^j^ = 0, implies 



Eliminating the /j,, v — components in Eq. fl57|) we can express the photon 
exchange in terms of its transverse and longitudinal polarizations, 



T = ^- 2 E (-l) A+ W(A)e(A)-j fc , (60) 

A=±1,0 



where 



£ (A = ±1) = --L(l,±i,0) , 

e(A = 0) = ^(0,0,1) . (61) 

We define the azimuthal angular orientation by fixing the hadronic vertex at 
tp = 0, while 

£' = (is(l-y)/yJ ± cos<pJ ± sm<p,£' z ) , (62) 
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where v = y£°. 

The following expressions will be given in the scaling limit, 

j ^ °° with x = — — — fixed, (63) 
[ v — > 00 2m N v 

where is the nucleon mass. We have then 



£ ± = £> ± = Q^l-y/y , (64) 

t = v(l-y)/y . 

The spin-averaged square of the scattering amplitude T in Eq. ( |6"0| ) can be 
written 



l^| 2 =^E(-l) A+A '^A^AA' , (65) 



A, A' 



where 



Explicitly, 



£AA' = i£j,-£*(A)j,*-s(A') = Ll, A . (66) 

spins 



L ±1>±1 = e 2 Q 2 [1 + (1 - /y 2 , (67) 
Afci*i = 2e 2 Q 2 e^(l-y)/y 2 , 
L ±1)0 = -V2<?Q 2 e™Jl^y'(2-y)/y 2 , 



£0,0 = 4e 2 Q 2 (l-y)/y 2 . 
The hadronic tensor is 

Hx>x = \j2r h - e*(A0 j fc • e(A) = F A * A , = tf_ A ,,_ A , (68) 

spins 

where the last equality applies if we choose H to be real, which is always possible 
by a choice of frame. With the expressions (|68|) for L\y we find for |T| 2 in Eq. 

GD, 

\T\ 2 = ^- 2 {^[l + (l-y) 2 }H u + (l-y)H 00 (69) 

+ (2 - y) v /2(l - y) cos<^ Re(H w ) + (1 - y) cos(2y?) i^} , 

which is the general expression for the hadronic tensor assuming only single pho- 
ton exchange between the lepton and hadron vertices (see for example Ref. ||17|| ). 



A2. Virtual Photon Wave Function. 
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The standard covariant expression for the 7 — > qq vertex at lowest order (Fig. 
A2) is 

iT^ M =u Xl (p 1 )Hee q 5 AB (-e(\)-j)]v X2 (p 2 )2n5(E y -E 1 -E 2 ) , (70) 

where E y = v is the (virtual) photon energy, J5j = \Jj> 2 + m? are the quark 
energies, A, B the quark colors, and we suppressed the 3-momentum conserving 
5-functions. The photon polarization vectors e(A) are given in by Eq. (|6l]) , and 
include the longitudinal (A = 0) polarization for virtual photons with q 2 = —Q 2 ^ 
0. 



, u , , ^ ^ q(p;A, ls A) 

q(p 2 ;A 2 ,A) 



Figure 6: Lowest order diagram describing the qq Fock state of the photon. 

We shall define the photon wave function V x XlM 5 AB by Eq. fl70|) with the 
replacement 

,0 

2vr5(£ 7 - - E 2 ) -> / dtexp[-it(E 7 -Ei -E 2 + ie)] , (71) 

implying 

^ AlAa <W = - ~ ee ^ , %( Pl )6(A) ■ 7^2) • (72) 

Physically, V^ A AlA2 is the lowest order amplitude for finding a qq pair at a given 
time (£ = 0), given a bare photon at t — —00. The time t of the transition 7 — » gg 
is summed over the interval — 00 < t < 0. 

The concept of a photon wave function is natural in the scaling limit fl63|) , 
where the mass Q of the photon is much less than its energy = v. In this 
limit, Eq. ( [72]) simplifies to 

^ AlA2 (z,p±) = -iee g ^ ' "^ (^(A) .7^) , (73) 

where z is the fraction of the longitudinal momentum carried by the quark, 

Pi = ^q + p ± , (74) 
P2 = (l-z)q-p± , 
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while 

e 2 = m 2 + z(l - z)Q 2 . (75) 
We shall use the helicity basis for the Dirac spinors, 

*m-#5;(*o W )' (76) 

where the 2-component spinors X\(p) are given in terms of the polar and az- 
imuthal angles 9, if of p as 



X tip) 



e"? sin(0/2) 

-e"^ sin(0/2) 
cos(#/2) 



In the scaling limit (^3]) we have 



. . vJ2z(l — z) 

x { P± e iA ^ Al ,_ A2 [z(A-2A 2 ) + (l-^)(A + 2A 2 )]-2m5 Ai2Al 5 AliA2 } , 

and 

^(z,px) = -iee q Al/Q[z $ 1 ]3/2 2A^ Al ,_ A2 . (79) 

It is useful to express the photon wave function in impact parameter space, 
through the Fourier transform 

V x XlM (z,r ± ) = /fe e ir -^ V x x ^(z,p ± ) . (80) 

The integral in Eq. (^) can be explicitly performed, for the photon wave functions 
in Eqs. (|78|) and ( [79]) , as well as for the other amplitudes discussed in the text, 
by using the identity 

oo y n+l 

du^—- I J n {u)=a n K n {a) , (81) 

where J n and K n are Bessel functions, and the identity applies for n < 3/2. The 
results can be written as 



V^±i(z,rj_) = — 1 ee q vJ2z(l - z) {-2m <5 Ai2Al (^Aa-KoO^jJ 

+ i e iA ^'5 Al ,_ A2 [z(\ - 2A 2 ) + (1 - z)(X + 2A 2 )] eK^er^)} 

V^ 2 {z,v L ) = ^ee q vQA[z{l-z)f /2 2\ l 5 Xl ^K {er A _) , (82) 
Zty 



-l 
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where ip r is the azimuthal angle of rj_. 

As shown in Section 2, the scattering amplitudes with one and two gluon 
exchanges off a massive target are, in the small x limit, proportional to the 
photon wave functions in Eq. (jS2p. 
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